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ABSTRACT 

The Tracy Puraping Plant consists of six large 
units with each pair of pumps connected through a wye branch 
to a 15 ft, diameter, reinforced concrete pipe which extends 
approximately 1 mile to the upper canal. Each pump is driven 
by 22,500-hp motor and delivers water at the rate of 767 cu-f1 
per sec. to the canal under a rate head of 197 ft. On the 
discharge side of each pump there is a butterfly valve which 
closes at a variable rate under the action of servomotor v/hen i 
power tripout occurs at the pump. 

A comprehensive study of the hydraulic transients 
for parallel' pump system for the Tracy Pumping Plant and for 
different modifications of the original plant involving a totei 
of 14 different cases has been carried out. For this a comput 
program has been coded in FORTRAN-IV and run on Dec- 1090 
system at I.I.T. Kanpur, The system has been analysed v/ithout ; 
a valve, with a butterfly valve at the pump, a non-return valv: 
at the pump, a combination of a non-return valve at the pump 
and an intermediate non-return valve at mid section of the 
discharge line. It has also been analysed with some suction I 

' ■ ' i 

length, increasing the length of discharge line, and doubling I 
the moment of inertia. The system has also been analysed with| 
an airchamber with a check valve upstream of the chamber, by j 



chsnQin^ "fch© diff ©xsn't psxsins'tsrs involvscJ in "tihs cornput©— • 
tions. All the cases, being analysed under this thesis, for 
the transient problems, are due to sudden power tripout* 

It has been found that maximum head is 1 ,4 times 
the rated head when, system has been analysed without valve. 

The system, analysed with butterfly valve has sho\/vn that 9 
percent reducation in maximum head is achieved. The transients 
for non-return valve at the pump are so severe and quick that 
special precautions will have to be taken. In the case of 
^^si stance dominated system, the maximum head rise is equal to 
steady-state rated head, and column separation may occur with 
subsequent rejoining and development of high heads. In case 
of gravity loading, doubling of moment of inertia has decreased 
the maximum head and increased the minimum head, and so increase 
in moment of inertia will be benificial. The system studied 
with air chamber has indicated that important design parameter 
for the design of airchamber is , the values of upsurge 

and downsurge decrease with increase in this constant. The 
0 f 1C .0 diame t er does not seem to have appreciable effect on 
upsurge and downsurge. 

The findings of the present study provide valuable 
suggestions for effective and rational design of surge controll- 
ing equipments. 



CHAPTER 1 


INTRODUCTION AND REVIEW OF LITExRATURE 
1,1 Introduction; 

Vfeter-hammer problems relating to turbo- 
machinery are of great practical importance. If the possible 
operating conditions of hydraulic-turbine and centrifugal- 
pump installations are compared, it soon becomes apparent 
that the pumps are subject to much wider and more involved 
variations than are the turbines, especially during the 
transient state of starting, -stopping, or emergency operation. 
In turbines the direction of flow and the direction of rotation 
are always the same, even in case of a breakdown of machine 
itself or trouble in the penstock and auxiliary equipment. 

Thus the machine performance always lies in the quadrant of 
normal turbine operation, and since its hydraulic characteri- 
stics are very well-Icnovm in this quadrant, it is completely 
straight-forward to predict the complete performance during 
any possible transient condition. On the other hand , under 
similar conditions with a pump installation, the flow can 
completely reverse direction, as can the rotation. The mach- 
ine in this case ceases to be a pump, and after passing through 
a zone of energy dissipation , becomes a runaway turbine. 

This great variation in performance gives rise to raany ques- 
tions, such as the runaway speed of the machine as a turbine-, 
the time of reversal, the magnitude of the accelerating 
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forces, the effect on the surge cycle in the discharge line 
maximum and terminal reverse rates of flov;, and so on. 
Unfortunately these questions are very difficult to answer, 
because, although the hydraulic performance of the machine 
is well known as long as it is acting as a pump, comparati- 
vely little study has aver been made of its performance as 
an energy dissipator or as a turbine. 

ThereforOj to explore these little-known regions 
of performance and to attempt to use the resulting informa- 
tion to answer hydraulic transient problems have their own 
significance. 

As the size of hydraulic-power units has increa- 
sed, more attention has been given to the design and location 
of auxiliary equipments for the modern hydroelectric plants 
such as governors, valves, surgetanks, and airchamber. The 
characteristics of a pipe system during transients depend 
upon the details of the pipe system and the cases of transi- 
ents. 

Butterfly valve covers a wide range of enginee- 
ring conditions and is suited to many hydraulic installations. 
There is reason to believe that this type of valve, because 
of its simplicity of design, its reasonable cost, and its 
inherent flexibility, has by no means reached the limit of 
its possibilities. The widespread use of butterfly valves , 
for controlling waterf low in large pipe lines proves their 
fitness and adaptability to this type of service. 
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V'/ater pumping rising mains usually have one or 
more non-return valves located along the main in addition 
to the non-return valve near the pumps. These valves may be 
referred to as intermediate non-return valves to distinguish 
them from the valve near the pumps. In many systems, inter- 
mediate non-return valves are provided with a feeling that 
it may help in controlling surge pressures, which may not 
be true. In pumping mains with surge protection devices 
such as airchamber, the effect of intermediate non-return 
valves on the transient pressures is of particular interest. 

Two extreme cases of pump loading are considered. 
First a relatively short line with almost all of the head used 
to lift liquid to a reservoir and second, for the same 
steady-state head and discharge, a very long pipe line in 
which substantially all head is used to overcome pipe line 
friction. Effect of moment of inertia of rotating parts is 
studied in two cases. Analysis of system with airchamber 
has also been carried out under different conditions. 

1.2 Review of Literature: 

A good amount of v\/ork has been carried out on 
hydraulic transient problems with the considerations of 
different auxiliary equipments. John Parmakian (5) studied 
the hydraulic transients with butterfly valve at the Tracy Puni« 
)ing Plant, According to him butterfly valves with a slow time 
of closure can be of no value in the control of transient 
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hydraulic conditions because the maximum rise, the maximum 
drop, and the reversal of flow in the pump discharge line 
occur within a short period. He also conducted tests inclu- 
ding following conditions, to check the adequacy of the 
pumps, discharge lines, and pressure control-equipment, and 
also to check the accuracy of water-hammer analysis: 

1, Two pumps supplying water to the same discharge line 
when a power failure occurred at both pump motors. 

2, Two pumps supplying water to the same discharge 
line when a power failure occurred at only one 
pump motor. 

3, One pump supplying water to the discharge line 
when a power failure occurred at the pump motor. 

Streeter (12) gave detailed information about 
different methods of analysing water-hammer problems, and 
he included the effects of type of pump loading on transient 
pressures, valve stroking applied to the pump pressures, etc, 

A number of studies (6) have been reported in 
literature on the effect of non-return valves and control 
valves (with set closure pattern and time ) on hydraulic 
transients. Murty, Shridharan and Prasad (6) studied the 
. effect of intermediate non-return valves on hydraulic tran- 
sients following power failure at two pumping systems and 
they illustrated the effects of (a) location of the 
intermediate non-return valve, (b) delay in the valve 



closure^ (c) bypass slze» (d) multiple intetmediate npn- 
feturn- valves, and (e) random delays in the closure of 
valves. They applied method of characteristics to Sit,V,i^ 
the effect and concluded that thete is no strong case for 
providing intermediate hon-return valve and where t^rovldedi 
they must be accompained by a bypass. 

Chaudhry (2) gave a widfi description of different 
boundary conditions for air-chamber and redticed the equations 
in such a form that it .can be used in computer analysis. 
Allievi (1) discussed the use of air chambers in pumping 
systems to control transients generated by |power failure to 
the pumps. Charts to determine the sire Pjf an .alrchamber 
for a pipeline to keep the maximum and minimum .puress^res 
within design limits are reported in the literature i(12-3). 
These charts may be used to determine the approximate size 
of chamber for a pipe line . 

1.3 Objectives of the Study: 

The major objectives of the study are : 

1) to develop a computer program for the transients in 
the pumping system due to power failure, 

ii) to study the effects of butterfly valves or non-return 
valves at the pump and intermediate non-return valves 
on the transients. 
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iii) to study transient in lift dominating versus 
resistance dominating pipe systems, and 

iv) to study the effect of an airchamber for mitigation 
of the transients and in this process to study the 
effect of various design variables relating to air 
chamber. 

1,4 Scope of the Study; 

Transient problems are broad and diverse, and 
include areas such as hydroelectric projects, pumped storage 
schemes, water~supply systems, nuclear power plants, oil 
pipelines, and industrial pumping system. In the present 
study transients in pumping systems have been studied. The 
transients in the Tracy Pumping Plant have been studied by 
Parmakian and Streeter (8). In the present study, a computer 
program has been developed for the analysis of transients in 
the pumping systems. With this program, the Tracy Pumping 
Plant has been analysed and the results have been compared 
with those obtained by Parmakian and Streeter. Subsequently, 
the effects of butterfly valve and non-return valve at the 
pump, intermediate non-return valve in the discharge line, 
comparative study of lift dominating versus resistance 
dominating systems, -'d^he effect of an air chamber to mitigate 
the transients, etc. have also been studied. 
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1.5 Significance of the Study s 

The work carried out will be useful in unders- 
tanding transients in, pump systems with different types of 
controls and safety measures, and in designing pumping 
systems to be safe under various possible transient 
conditions. 

1.6 Organisation of the Report; 

The study is reported in the following sequence; 

i) The analysis and comparison of results for Tracy 
Pumping Plant, description of sequence of events after 
power failure, derivation of differential equations of 
rotating mass and boundary conditions are presented 

in Chapter 2, 

ii) The analysis of hydraulic transients with modifications 
in the Tracy Pumping Plant has been done in Chapter 3. 

iii) General boundary conditions , method of solution of 
transients and results for airchamber are presented in 
Chapter 4. 

iv) Summary, conclusions and. suggestions for future study 
are given in Chapter 5. 



CHAPTER 2 


ANALYSIS OF TRANSIENTS OF THE TRACY PLWING PLANT 

2.1 Sequence of Events Following Power Failure: 

In the analysis of piping systems which contain 
turbomachines it is necessary to take into account the chan- 
ges occurring in the characteristics of the units as their 
rotational speeds change. During a pump start-up, the 
discharge valve is usually kept closed to reduce the elect- 
tical load on the pump motor; and as the pump speed reaches 
the rated speed, the valve is gradually opened. It is the 
usual practice, in a normal pump~stopping procedure, that 
the discharge valve is first closed slowly, and then the 
power supply to the pump motor is switched off. In cases 
of large pipelines designed for low heads, a bypass to the 
suction reservoir may be used in starting axial flow units. 

In the above procedure of starting and stopping pump motor, 
the speed of the pump remains almost constant during the 
transients in the piping system. Transients caused by emer- 
gency pump operations are usually severe, and the pipelines 
should be designed to withstand positive and negative pressures 
caused by these operations. There are many operating condi- 
tions at large motor-driven centrifugal pump installations 
which are capable of producing substantial pressure changes 
in the discharge line. Of these conditions , one of the most 
important is rapid deceleration of the pump motors because 
of a power failure. 



9 


When power to a pump fails while it is lifting 
liquid to a reservoir, the follovvdng events take place in 
the absence of the closing of a discharge valve.. The pump 
speed reduces since pump inertia is usually small compared 
to that of the liquid in the discharge line. Because the 
flow and the pumping head at the pump are reduced, negative 
pressure wave propagates downstream in the discharge line, 
and positive pressure v;ave propogates upstream in the suction 
line. Flow in the discharge line reduces rapidly to zero and 
then reverses through the pump even though the latter may 
still be rotating in the normal direction. In this condi- 
tion (i.e,, v;hen there is reverse flow through the pump 
while it is rotating in the normal direction), the pump is 
said to be in the zone of ’’ energy dissipation Because 

of the reverse flow, the pump slows down rapidly, stops mom- 
entarily, and then reverses, i.e. the pump is now operating 
as a turbine. The pump speed increases in the reverse dire- 
ction until it reaches the runaway speed. With the increase 
in the reverse speed, the reverse flow through the pump is 
reduced due to choking effect, and positive and negative 
pressure waves are produced in the discharge and suction 
lines, respectively. 

Vifhen the load on the pumping system is primarily 
due to fluid friction, as in the case of a long discharge 
line, the flow is retarded very slowly and no reversal may 
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occur at all. The danger to the discharge line in this 
case is the development of column, separation from low/ 
pressures. If the pipeline profile is such that the trans- 
ient-state hydraulic grade line falls below pipeline at any 
point, vacuum pressure may occur, and the water column in 
the pipeline may separate at that point. Excessive pressure 
will be produced when the two columns later rejoin. During 
the design stages, the possibility of water-column separation 
should be investigated, and if necessary remedial measures 
should be taken, 

2,2 The Tracy Pumping Plant: 

The Tracy Pumping Plant is located in the central 
valley project near Tracy, Calif, USA. The general arrangement 
of the plant is shown in Fig, 2» The plant consists of six 
large pumping units with each pair of pumps connected through 
a way branch to a 15-ft. diameter, reinforced concrete pipe 
which extends approximately 1 mile to the upper canal. Each 
pump is driven by a 22,500-hp motor and delivers wfater at the 
rate of 767 cu-ft. per sec to the canal under a rated head of 
197 ft. On the discharge side of each pump there is a butterfly 
valve which closes at a variable rate under the action of ser“ 
vomotor when pov/er failure occurs at the pump. There are 
several operating conditions at the Tracy Pumping Plant that 
can cause substantial pressure changes in the discharge lines. 
The most serious of these operating conditions is the sudden 
stoppage of the motors as a result of a power failure. 
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The transients were analysed by Parrnakian and Streeter and 
the time history of the pressure, discharge and speed changes 
caused by a simultaneous power failure at the motors of two 
pumps on the same discharge line were presented (Fig. 2a). 

It can be seen that maximum head rise at the pump is 
1.18 Zone of normal operation is upto 9.8 sec, zone of 

energy dissipation is from 9.8 sec. to 16 sec., and 
zone of turbine operation is from, 16 sec. The results are 
shown upto 32 sec. The Tracy Pumping Plant analysed without 
butterfly valve gave maximum head ^rise as 275.8 feet as stated 
by Parrnakian. 

2.3 Perform.ance Characteristics of a Pump; 

The discharge of a centrifugal pump depends upon 
the rotational speed N, and the pumping head H; and the 
transient - state speed changes depend upon torque T, and 
the combined mom.ent of inertia of the pump, motor, and liquid 
entrained in the pump impeller. Thus, four variables- namely, 
Q, H,N, and T- have to be specified for the mathematical repre- 
sentation of the dyanamics of a pump. The curve showing the 
relationships between these variables are called the pump 
characteristics. There are several methods to store these 
characteristics for computer analysis. For storing pump 
characteristics in 'a digital computer, the method used by 
Marchal, Flesch, and Suter (12) appears to be most suita- 
ble and is used herein. 
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Although pump-characteristics data in the pumping 
zone are usually available, only very limited data are 
generally available for either zone of energy dissipation 
or the zone of turbine operation. If the complete characteri- 
stics data are not available, then the characteristics of a 
pump having about the same specific speed may be used as an 
approximation* Homologous theory permits results of model 
tests to be used; and by using dimensionless, homologous 
relationships one can store the characteristic curves in the 
computer in very few units of storage space. To utilize these 
characteristics, by interpolation for gate position or vane 
angle, one finds three points on the characteristic curve in 
the vicinity of the desired conditions. A parabola is fitted 
through these three points for calculation of boundary condi- 
tions at the particular instant under consideration. 


Data for prototype pump characteristics are 
obtained from model test results by using homologous relation- 
ships. Two pumps are considered homologous if they are geome- 
trically similar and the stream flow pattern through them is 
also similar. Then they will satisfy the conditions for 
dynamic similarity except for viscous effects. The homologous 
conditions are given by 


H 


nr 


= constant 


(Eg. 2.1) 


and 


= constant 
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In other v;ords, these ratios must be the same 
for any of the homologous series of units when they are 
operating in a dynamically similar manner. Vi/hen studying 
transient effects in a given turboraachine whose size is 
fixed and known, diameter of the impeller may be included 
in the constants; hence 


= constant 
— ^ = constant 


(Eq. 2.2) 


Now if the ratio of li, Q, and N to their rated 
values Nj^ are used to yield a dimension less presenta- 

tion let 

h = , V = , and a = ^ (Eq.2a) 

The dimensionless-homologous relations are 



constant , — = constant 

a 


(Eq.2b) 


Since a becomes zero while analysing transients 

for all four zones of operation — ^ becomes infinite. To 

h h . 

avoid this, the parameter — k — instead of may be 

(a'^+v'^) a 


used. 


The signs of v and a depends upon the zone of 
operation. In addition to the need to define a different 
characteristic curve for each zone of operation, ^ becomes 
infinite for v=0. To avoid this, a new variable 0 may be 


defined as 



(Eq. 2c) 


0 = tan“^ ^ 

and then characteristic curve may be plotted 
2 2 

between 0 and h/(a 4-v ) , By definition, © is always finite, 
and its value varies between 0° and 360° for the four zones 
of operation given in Table 2.1 below; 

TABLE 2.1 ZONES OF PUaIP OPERATION 


Zone of operation 

V 

a 

Range of © 

Pump 

4- 

+ 

0°<©<90° 

Energy dissipation 

- 

4- 

90° <0 <180° 

Turbine 


- 

180° <©<270° 

Turbine energy 
dissipation 

4- 

— 

270° <©<360° 


If T be the shaft torque required for given 

speed and discharge ratios, the dimensionless-homologous 

T 

relation for torque ratio, p = are 


^ = constant 

a 

and “ = constant 


(Eq.2d) 


Similar to the head characteristic curve, the 
torque characteristic curve may be plotted between 
p/(a^+v^) and 0. 
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2.4 Equations of Conditions Imposed by Pump: 

, The pump characteriestics curve can be used if 

we have discrete points on these curve at equal interval of 
0, between the range 0=0° and 0 = 360°, are stored in the 
computer. Each segment of these curves between points stored 
in the computer may be approximated by straight lines. If 
sufficient number of points are stored in the computer then 
error introduced by approximating these curves by segmental 
straight lines is negligible. 

Given that variables a,v, h and p are known at 
i-th time step it is needed to calculate these variables at 
the end of time step denoted respectively by oc ,v ,h , and p . 

)r )r 

Since small time steps may economically be used 
in carrying out the transient calculation on a large digital 
computer, the speed change of a turbomachine may be very 
closely estimated for a particular time step if one extra- 
polates the speed and discharge ratios to the mid point of 
a small time step, calculates the torque for these conditions, 
and then determines the speed change at the end of the time 
step from torque and moment of inertia of rotating members. 

To determine the values of these variables it is in turn 
necessary to determine the equation of the segment of pump 
characteristics corresponding to a and v , As a first esti- 
mate, their values are determined by extrapolation from the 
known values from the previous time step i.e., 
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a 


e 


+ 


i-1 


V. 




^^i-1 


(Eq.2.3) 


in which and are estimated values at the 

end of i— th time step, and refers to known values at 

the beginning of i-th time step, and£v.aj_^ and.£iV^__j^ are 

changes in these variables during the (i--l)th time step. 

Now, the grid points on either side of © = tan""^ (a^/v^) are 

searched, and the ordinates h/ (a +v ) and p/(a‘'+v"^) for 

these grid points are determined from the stored values. From 

these, the constants aj|^,a 2 »a 2»®4 the equation of the 

segmental straight line are determined. Now, assuming that 

the points corresponding to a ,v ,h and 6 lie on these 

p p . p p 

straight lines, then 


P P 


= a, 4-a^ tan~^ 


a. 




V. 


(Eq.2.4) 


and 




R. 


2 , 2 
“p'-^P 


*•1 

a-j+a^ tan~ — ^ 
3 4 Vp 


(Eq.2.5) 


in v/hich a 2 Ta 2 and a^ja^ are constants for the 
straight lines representing, respectively, the head and torque 
characteristics. 

Boundary conditions at the turboraachines may 
include combinations of units in parallel and in series 
cases in which some units lose power while others continue 
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to operate, and action of valving. Since the pipeline 
water-hammer is conveniently calculated by the method of 
characteristics, the inclusion of various appurtenances 
along the pipeline to control the transients may be included 
or omitted from the program as an option that in no way 
changes the more complex boundary conditions at the turbo- 
machines. 


Refering to Fig. 2.1 the following equations can 
be written for the total head at the pump; 

‘V “ “sue + V ”pv- (aq.2.6) 

i,l 

in v/hich, 

^suc ~ height of the liquid surface in the 
suction reservoir above datum, 

Hp = pumping head at the end of time step 

Hp^ = head loss in the discharge valve. 

The value of head loss is given by 


^Pv “ ^v^ ^P.. /I^P.. 

1,1 1,1 


( Eq. 2.7 ) 


v 


in v;hich = coefficient of head loss in the 
alve, and Qp. = discharge at section (i,l). 


2.4.1 Characteristic Equation for Discharge Pipe: 


As the suction line is short, it may be neglected 
in the analysis, The negative characteristic equation may 
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be written along AP, (Fig. 2.2) 



'' “a - 


llJlt 


2.D.A 


X 





^ X H„ 
a P 


or Op = Cn + Ca. Hp (Eq.2.8a) 

where Cn = Q^- x - 2, d/A ^ I ^A^ 

and Ca = 

a 

Reffering to Fig. 2.1 characteristic equation for 
section (i,l) may be written as, 

Qp. = Cn+Ca.Hp (Eg. 2.8) 

*^i , 1 ^i , 1 

2.4.2 Continuity Equations; 

Since there is no storage between the suction 
reservoir and section (i,l) it can be written as 

Qp_ = Qp (Eq.2.9) 

i , 1 

in which Qp = flow through the pump at the end of 

time step. 

2.5 Differential Equation of Rotating //lasses and Calculation 
of Speed Change; 

When the power to, the pump motor is interrupted 
or suddenly cut-off, the deceleration of the pump at any 
instant depends on the flywheel effect of the rotating |>arts 



Hjr€r§iifh tr9i§, //## 
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of the pump and motor and the instantaneous torque exerted 
by pump impeller. For a rotating system, the accelerating 
torque is equal to the product of the mass moment of inertia 
of the rotating system and the angular acceleration. Following 
a power failure at the pump motor, the decelerating torque of 
the rotating' system corresponds to the pump-input torque prior 
to power failure. If the decelerating torque is taken as 
positive, 

T -I dt g dt (iiq,2.10) 

in which 

T = the pump input torque corresponding to a given 
speed and head, in pound-feet 
I = the mass moment of inertia of the ratating 
parts, in pound-feet squared 
0 ) = the angular velocity of pump and motor shaft, 
in radians per second 
t = the time in seconds 

A = the flywheel effect, or moment of inertia in 
terms of the weight of the rotating parts of 
the motor, pump and entrained water, in pound- 
feet square 

g = the gravitational acceleration, in feet per 
second squared. 

Because the pump torque varies with the speed, 
Eq.2.10 is applicable for specific instants of time only. 
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For a small time interval (At = Eq.2.l0 can be 

approximated as 



>L _ 2 tz a ^ 1’“^'^2 

g at ” 60.g ^ At 


(Eq.2.Jl) 


in which N is the pump speed, in revolutions per 
minute. Eq.2.11 can be written as 

IS.gTj, 

“l‘‘°^2 “ u'a Nq (x^q.2,12) 

in which a = N/Nj^, p=T/Tj^ and the subscript R is 
used to denote the rated conditions. ;declerating torque 
at the rated head and speed of the pump is 


T = ^O'^HrQr 

2.n.Nj^Tip^ 


(Eq.2.13) 


in which y is the specific weight of water, in 
pounds per cubic ft. ; H represent the pressure head for 
surge conditions measured above the pump-intake; watersurface 
elevation, in feet; Q is the pump discharge, in cubic foct per 
second and x] denotes the pump efficiency. Then 


- “2 = 

450.g,Y,HpQn 

in which K,= £LJi = 

^ 9 9 

. A T}j^ Nj 


91,600 

A 


Eq.2.14 can be rewritten as 
ttp - C^. pp = oc + C^» P 

-15 Tp. A t 


in which = 


— ^ 

•n.wR N, 


'R 


( £:q .2. 14 ) 
(Hq.2.15) 


(Eqn.2.16) 


(Eq.2.17) 
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2.6 Boundary Conditions for Parallel Pumps~ Short Suction Line: 

The continuity equation for parrallel pumps is 

Qp ~ ^p*Qp 

"^i,l ^ ^ 

in which np = number of parrallel pumps. 

Depending upon the length of the suction line, 
boundary conditions for parallel pumps may be divided into 
the following two cases: 

i) short suction line 
ii ) long suction line 

In case of short suction line, the water-hammier 
waves are neglected. By eliminating Hp , Hp^ and Qp,^ from 
eqns. 2. 6, 2. 7, 2. 3, and 2.9 and by using and as reference 
values, the resulting equation may be written as 

Qji^.Vp = Cn Ca.Hj^.hp -- Ca.C^.Q^ ^pl^p! (cq*2»i‘9’} 

Now we have four eqns.“i.e. ’-^q. 2.4,2.5,2.26, and 2.19 in four 
unknowns ap,hp,Vp,Pp. 

By substituting for hp from eqns. 2.4 into 2,19- a:xi 
for Pp from eqns. 2, 5 .into eqn.2.16 and simplyfying, we 
have 

= Ca • Hj^. (cXp +Vp ) + Ca.Hp a 2 (ctp + Vp )tan 

- C^.Vp. np - Ca.C^. C? Vp|vp| +Cn+Ca.H^y^= 0 {«j.2.20). 

• R 
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1 

b 

II 

CM 

- / 2 2 X ^ / 2 2x 

, -1 

tan —=■ -a 

-c^.p =0 




(Hq.2.21) 

. 57; = 

Cg • Hjp^ ( 2 • 3 • v^*" 32 ♦ ^ 'tisn 

T a 

"p- 

% 


— 2. Ca . 0^ . 1 Vp 1 


(Eq.2.22) 

aPi _ 

■bap 

Ca.Hp(2apap+a2*Vp+2. a2»a:p tan“^ 

CTt 

(iq.2.23) 

1!2_ 

dap 

1- C^( 2.a3.ap +a4,Vp+2. a^. ftp tan — 

{Eq.2.24) 

and 

■^p 

C^(~2.a2»Vp+ a^.ap-2,a^.v^ tan“ 

1 «p>. 

^p 

(Eq.2.25) 


Let and be the initially estimated values 

of solution, which may be taken equal to a. and v, as determi- 

G 0 


ned from Eq. 

Then a 

better solution of Eqs.2.20 

and 2.21 is 

“P 

= 4^> 

+ 5ap/2 

(Eq.2,26) 

(2) 

^P 

II 

< 

H* 

+ 5vp/2 

(2q*2.27) 


subscript (l) indicates estimated values and subscript (2) 
indicates valves after first iteration. 

in which 


AF AF 



(Eq.:>.2=8i 







2* da. 


- F 


and 


s; 


dF. 

1* aZ 


V, 


dFj, 


aF^ dF^ dF^ 


(Eq. 24 i 29 ) 


da. 


da,^* dv, 


Lp WW.p UVp 

2,7 Procedure for Solving the Transients: 

The general procedure for solving the transient 
problem is to first set up the steady-state case for the given 
elevations of reservoirs, and pump and pipeline characteristics, 
computing steady-state head steady-state discharge Q^, and 
the velocity, and elevation of hydraulic grade line at equally- 
spaced sections along the suction and discharge pipes. The 
basic data required for analysis of the problem are of Tracy 
Pumping Plant given in Table 2,2 and the head and torque data 
required are given in Table 2,3, The tabular data of Table 2,3 


are (a/v or v/a) vs ( 


h 


3 or “ 
a~ V" 

To use Table 2,3 in the computer it is converted to Table .2.4 


and ( ~ or vs (^ 2 °^ 

^ ^ “ a V 


in which ( 7 or 7 ) vs ( and ( 7 or 7 ) vs ( 

values are stored. By this computation , theta (© = tan“^a/v) 
values were found to be at unequal intervals. So this table 
is further changed to Table 2,5 to obtain tho values of © at 
equal intervals . After obtaining Table 2,5 , the speed change 
is calculated from Eq.2.3 Section 2.4, 


A sketch of sectional elevation and plan of pump 
and discharge line is given in Fig. 2. It is clear from the 
figure that the upstream boundary condition is imposed by 
pump characteristics and dounstream boundary condition is 
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Basie Bata and Results for different eases without Air Chamber 
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Table 2.4 Converted head and torque data for pump of 
specific speed 1,800 Cgprn units) 
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Table 2.5 Pump Characteristics Data at Equal Interval of © * 


©= tan"^ 

a 

V 

NS = 1800 

(gpm units) 



h 

B 

9 

h 

- 3 


"'7^2 
a +v 

2 2 
a +v 

a. +v 

2^ 2 
a +v 


-0.556 

-0.372 

140 

0.957 

0.598 

5 

-0.478 

-0.299 

145 

0.920 

0.671 

10 

-0.386 

-0.209 

150 

0.886 

0.740 

15 

-0.284 

-0.089 

155 

0.860 

0.798 

20 

-0.181 

0.026 

160 

0.830 

0.839 

25 

-0.090 

0.137 

165 

0.799 

0.865 

30 

-0.076 

0.309 

170 

0.766 

0.881 

35 

-0.186 

0.324 

175 

0.730 

0.883 

40 

0.356 

0.416 

180 

0.692 

0.865 

45 

0.500 

0.500 

185 

0.655 

0.828 

50 

0.643 

0.555 

190 

0.617 

0.782 

55 

0.783 

0.595 

195 

0.589 

0,745 

60 

0.909 

0.614 

200 

0.563 

0.678 

65 

1,020 

0.615 

205 

0.538 

0.589 

70 

1.113 

0.600 

210 

0.519 

0.496 

75 

1.188 

0.570 

215 

0.511 

0.403 

80 

1.244 

0.535 

220 

0.507 

0.316 

85 

1.276 

0.494 

225 

0.506 

0.228 

90 

1.288 

0.450 

, 230 

0,520 

0.143 

95 

1.279 

0.397 

235 

0.541 

0.074 

100 

1.266 

0.365 

240 

0.567 

0.006 

105 

1.244 

0.344 

245 

0.593 

-0. 050 

110 

1.207 

0.331 

250 

0.615 

-0.118 

115 

1,157 

0.333 

255 

0.632 

-0.224 

120 

1.112 

0.351 

260 

0.642 

-0.360 

125 

1.074 

0.391 

265 

0.646 

-0.684 

130 

135 

1,035 

0.996 

0.444 

0.520 

270 

0.634 

-0.684 
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due to constant head reservoir. In case 1. of the study 
there is no valve just downstream of the pump and the 
transient is due to failure of two pump motors. 

The suction pipe length is zero in this case and 
length of discharge line is 5130 ft. The plant consists of 
six large pumps with each pair of pumps connected through a 
wye branch to a 15-ft. diameter, reinforced concrete pipe 
which extends approximately 1 mile to the upper canal. 

Each pump is driven by a 22, 500~hp motor and delivers 
water at the rate of 767 cu-ft, per sec. to the canal under 
a rated head of 197 ft. The transient due to simultaneous 
falure of both the pumps are studied. 

At the upstream boundary we have four Eqs, 2,4, 2. 5-, 

2.16 and 2,il9 in four unknowns ap,Vp,hp and Pp . These four 
equations are reduced to two non-linear equations in two unkno- 
w/ns, ap and Vp. These two Bqs.2.20 and 2.21 are then solved by 
Newton-Raphson method , in which a solution of the e-quation is 
first guessed, which is then refined to a required degree of 
accuracy by successive iterations, I-Iaving determined ap and 
Vp, it is verified whether the segment of the pump characteri- 
stic used in the computation corresponds to ap and Vp. If it 
does not, then a and v„ are assumed eaual to a^ and v^. if 
the correct segment was used, then hp and Pp are determined 
from 2qs, 2.4 and 2.5, Hp and Qp are determined from the 

equations Hp=hpxHp and Qp=VpXQp; and Hp and Qp from 

r , 1 X , 1 
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Eqs. 2.6 and 2.9 . The values of a and v are initialised for 
the next time step ( i.e. a=ap and p=Pp), and the solution 
progresses to the next time step. To avoid an unlimited 
number of interations in the case of divergence of solution,’ 
a counter m.ay be used so that the computations are stopped 
if the number of interations exceed a specified value. 


For computations at interior points, let conditions 
at time t= t^ be known, these are initial steady~state condi- 
tions at time The unknown conditions at t^+ Dt are to be 

determined. Referring to Fig. 2.3 , positive characteristic 
equation may be vi/ritten AP, 


Qp— 


q. A u f . Dt r. I r, I 


a P 


or Qp= CP-Ca.Hp 
where, CP= Q^+ 


(Eq.2.30) 


f . Dt 

• ^A *“ 27071 ""A '"‘A' 


' Qa IQaI ■ 


and Ca 




Similarly negative characteristic equation may be 
written along the characteristic line BP, 


«P = QbI°bI + "r • Hp 




OX Qp *— Gn4"C3 1 Hp 


(Eqs. 2. 31) 


where, Cn - Qg - 


a 



f. Dt 
2.D.A 




and 
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In Fig. 2.3 characteristic lines AP and 3P have 
slopes + l/a. Physically, these lines represent path trans- 
versed by a disturbance. In 3qs. 2. 3Gnnd 2.31, there are two 
unknowns, namely, Hp and Qp. The values of these unknows can 
be obtained by simultaneously solving these equations, i.e», 

Qp = 0.5 (CPi +Cn) (Eq,2*32) 

Knowing Qp,the value of Hp can be obtained either by 
Eq#,2.3C or 2.31. Thus using Eqns. 2.30 and 2.31 conditions 
at all interior points at the end of time step can be deter- 
mined. 

For the calculation at downstream boundary, Hp 
is known and it is equal to downstream reservoir level, 

HRES. Qp can be determined using positive characteristic 
Eq. 2.30. The flow chart of Fig, 2b, and Fig.2cillustrates the 
procedure for calculation, 

2.8 A Computer Program for Transient Analysis of Pumps in 
FOi 1 TRAN Language: 

A computer program for analysis of the transients 
in a pipeline caused by pumps is available (2). It is suita- 
ble for analysis of easel Tracy Pumping ' Plant without valve. 
Hence above program was implemented on DEC-1090 system of 
I.I.T. Kanpur and the implementation was tested and verified 
with data given by €haudhry(2). The results agreed with 
the program output given in the reference and hence the imple- 
mentation of the program is considered satisfactory. 


1C : IPRiNT 


Flow CFocrt -^OY Ol SeYLes S^s-fc^vn 














FLoW CKcXyt OUJT^ cl COT^ci cfciori5 Foy RjCyY\jQ 
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2*9 Results for Case 1: 

The results obtained 

Using the data for the Tracy Pumping Plant without 

and 

valveZ.analysed with the computer program of Chaudhry (2), The 

results obtained are plotted as, (h vs t), (v vs t), and 

the 

(a vs t ) in Fig. 2,4. Fig, 2. 5 represents, /^transient pressure 
gradient in the discharge line as a result of power failure. 

It is found from the graph that maximum head obtained ia l,4Hj^, 
which is the same as obtained by Parmakian (6). 

2,10 Results of the Present Study and Comparison of the 
Results Presented by Streeter: 

Case 2, has the same plan, sectional elevation and 
basic data as for case 1 except for the fact that there is a 
butterfly valve just downstream of the pump. The method of 
solution is same as in case 1. except after 14 sec. Butterfly 
valve closure curve represented in Fig, 2, 6 has to be taken 
into account. The values ofC for different time steps has 
been interpolated by using parabolic interpolation. 

The graphs plotted for this case are, (h vs t), 

(v vs t) and (a vs t ) as shown in Fig. 2.6, Another graph 
has been plotted for this case as shown in Fig. 2,7. repre~ 
senting transient pressures in the discharge line as a result 
of power failure. It can be seen from Fig. 2.6 that maximum 
head is 1.3 This value indicates that there is decrease 

of of the maximum head as obtained without butterfly valve. 
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According to Parmakian (9) this decrease may be up to 20y< 
if proper closing of butterfly valve is done. The pressure 
drop at the pump is nearly the same as that Vi/hich would result 
without valve and will occur in 3,5 sec. after power failure 
of the pump motors; a complete reversal of flow in the disch- 
arge line will follow in approximately 9.8 sec. Therefore , 
valves with a slow time of closure can be of no value in 
controlling hydraulic transient conditions, because the 
maximum pressure rise, the maximum drop in head, and the 
reversal of flow in the pump discharge line, will occur within 
a short period of time. It is found that difference between 
computed and observed results in normal zone of operation is 
negligible for which correct data has been supplied by the 
manufacturer. Since assumed data has been used in energy 
dissipation and turbine operation zone, difference is appre- 
ciable. But plotted result for (v vs t) graph obtained in 
this study is much closer to that obtained by Parmakian (9). 
Comparing Fig. 2.5 and Fig. 2, 7 it is observed that region 
of negative pressure has got reduced. 
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FAILURE FOR CASE 2 


CHAPTER 3 


ANALYSIS OF TRANSIENTS OF THE TRACY PLi/'iPING PLANT VrtTH 
JVDDIFICATIONS IN THE SYSTEM 

In order to study the transients at Tracy Pumping 
Plant with modifications in the system the following special 
cases were studied in detail viz,, case 3 a non-return valve 
at the pump, case 4 a non-return valve at the pump and 
intermediate non-return valve at the middle of discharge 
line, case 5 lift dominating with two different moments of 
inertia of rotating parts and case 6 resistance dominating 
with two different moments of inertia of rotating parts. The 
detailed data for these cases are also shovm in Table 2.2 

3.1 Case 3 Tracy Pumping Plant with a non-return valve: 

Case 3 has the same plan, sectional elevation and 
basic data as for case 1 except for the fact that there is 
a non-return valve just downstream of the pump. The method 
of solution is same as in case 1, in addition boundary condi- 
tion at the pump changes when flow reverses. As flow reverses, 
the check valve gets closed and pump characteristics vanish 
from the computations. 

The graphs plotted for this case are, (h vs t), 

(v vs t), and (a vs t) as shown in Fig. 3.1, which represent 
transient as a result of power failure. Another graph 'has 
been plotted for this case shown in Fig. 3.2, which shows the 



LS = 0-0 LD = 5130 (FEET ) 

WR2 = 4040000 00 (LB-FEET2) 
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transient pressure in the discharge line as a result of power 
failure. It is found from Fig, 3,1 that head rise at the 
pump has gone as high as 1.9 Hp,, and minimum pressure is 0.14Hp. 
As discharge line is not too long and frictional effect is less, 
it is a case of gravity loading. Thus, for the gravity loading 
with all pumps failing, with undamped non''-return valve, very 
severe fluctuating transients occur as seen from Fig. 3.1, the 
water in the discharge line acts as a liquid spring with the 
hydraulic grade line at the valve oscillating about downstream 
reservoir* Since transient are so severe for complete failure 
of pump motors for gravity loadings, special precautions should 
be taken. Fig. 3.2 shows that the region of negative pressure 
for this case has got reduced as compared with case 1 and 
case 2, so chances of column separation has also got reduced, 

3,2 Case 4 Tracy Pumping Plant with two non-return Valve; 

The problem studied under case 4, is exactly same as 
in case 3, but in addition to case 3, there is an intermediate 
non-return valve at the mid section of the discharge line. For 
analysis of the problem the discharge line has been treated as 
two pipes in series, the length of the first pipe is considered 
upto mid section and second pipe starts from the mid section 
to the downstream reservoir. In the analysis of the problem 
three sections are of particular interest, first at the pump, 
second the section just before intermediate non-return valve. 

The boundary condition at the mid section' is changed when 
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TWO PUMP MOTORS 


47 



CC 

3 UJ 

to to 
to < 


i33d N1 dWnd dO 1 ) 3A0aV iH913H 



48 


flow reversed, as intermediate non— return valve gets closed 
at that instant. The positive characteristic equation has 
to be used for section second and negative characteristics 
equation has to be used for section three. The solution 
procedure is same as for the case 3. 

The variation of head at these sections are plotted 
as shown in Fig. 3,3, and it is observed that the fluctuations 
in head are not so severe and sharp as observed in case 3. 

The maximum pressure developed at the pump is nearly rated 
pumping head. At the second section the maximum head develo- 
ped is the steady-state head, but head rise at the third 
section is as high as 1,39 H^,, nearly the same maximum head as 
in case 1. Special auxiliary equipments has to bo provided 
whenever intermediate non-return valve has to be used. Another 
graph has been plotted as shown in Fig. 3,4, which shows the 
transient pressure, in the discharge line as a result of power 
failure. It is observed that region of negative pressure has 
been increased as compared to other cases. 

3.3 Pumps with Long Suction Length; 

If the suction line is not short compared to 
discharge line, then water-hammer in the former has to be 
considered in the analysis. Therefore, characteristic equa- 
tion for the suction line has to be included as follows 
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Refering to Fig, 3.5 


Hp 

= H, - H 

i+1,1 i , n+1 

(Eq.3.1) 

Op 

= CP-Ca4.H.p 

^ i,n+l 

(Eq.3.2) 

i , n-j-1 

^i+1,1 

= Cn+Ca . . H p 

^i+1,1 

{Eq.3.3) 

Clp 

^i , n+1 

= Qp = np, Qp 

^i+1,1 ^ ^ 

(Eq.3,4) 


multiplying Eq,3,2 by and Eq.3.3 by Caj^^ and substituting 

f03r Qjp and Qp ,and adding 

i +1 y 1 i ^ n+1 

= Cn.Ca^+CP.Caj^^j^.Hp (Eq.3.5) 

by using and Hp as reference values, Eq.3.5 may be written as 
^ ^ ^ Qp^* Vp""Cn • Cajj^*“CP» (Eq.3.5) 

Ca^.Ca^_^jp.Hp 


by eliminating hp from Eqs, 4 and 3.6 

= a^(Q:^ +v^) + a2(otp +v^) tan”^ ^ - C^.Vp+Cg=0 




(Eq.3.7) 

in which Cy 

np.Qp(Caj^+Cai_j^l) 

Ca^* Ca^_i_j_.Hp 

Cn# C 3 jj^’f'CP • 

Ca^.Ca^^p.Hj^ 

(Eq.3.8) 


and Cg 

(Eq.3.9) 

dFp 

"'San 

—1 '^P 

~ 2 . * oCp+2 . a^* ^p • tan ^ * ^F 

(Eq.3,10) 


-1 “P 

5 ^ = 2.a,.vp + a.a^.Vp tan - - a^.cc^-C^ 



and for this case 
2.24 and 2.25. 



and 




are given by 4q$,2.2i> 


In the study of transients in pipelines due to 
failure of pumps two different cases can be identified, one 
in which the resistance of the pipeline is small and hence 
the lift J-S dominating, and the other where the lift is 
small and resistance is dominating since the transients in 
the pipelines under these two cases are very different it 
is proposed to study them separately, 

3,4 Case 5 Pipes with Lift Dominating; 

The system studied has 1282.5 ft. of suction 
pipe and 3847.5 ft. of discharge pipe with a head loss due 
to resistance under steady~state of 3,9 ft. which is small 
in comparison to the gravity lift of 197.1 ft. 

3.4,1 Case 5 a. Small Moment of Inertia of Rotating Parts: 


This case takes into account the length of suction 
pipe, and as a result characteristic equation for the suction 
line has to be taken into account. 


expression for F^, , 


As suction length has not been neglected, the 

taken for case 1 is now 

dF, 


and 


dF 

dv, 


changed. The expressions for these variables 

dF P 

and ^ are defined by the Eqs. 3.7, 3.10 and 3.11 respectively, 

These expressions has two constants Cj and C 0 which have been 
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defined by the Eqs. 3.8 and 3.9, The solution procedure 
adopted is same as for case 1. 

The graphs plotted for this case are> (h vs t)^ 

(v vs t), and (a vs t) as shown in Fig, 3,6, which represent 
transient condition as a result of power failure. The maxi'" 
mum steady-state head is 2Q1 ft. which goes as . high as 

275.2 ft. 

3.4.2 Case 5b Large Moment of Inertia of Rotating Parts 5 

This modification in the Tracy Pumping Plant has 
been' studied under case 5b. of the computer program. This 
case is exactly same as case 5a, except moment of inertia 
of rotating parts has been doubled, which is being equal to 
8080000 pound feet squared. The equations involved are same 
as for case 5a and solution procedure adopted is same as for 
case 5a. 

The graphs plotted for this case are, (h vs t), 

(v vs t), and (a vs t ) as shown in Fig, 3,7, which represent 
transient conditions as a result of power failure. The 
maximum steady-state head is 201 ft. which goes as high as 
251.25 ft. The effect of moment of inertia can be seen from 
Fig. 3.8. It is found that doubling the moment of inertia 
lowers the maximum head developed along the discharge line 
but increases minimum head along the discharge line. The 
friction drop is very small for case 5a and hence 
loading is due to gravity lift. It is observed that increase 






NO VALVE 

LS = 1282-5 FEET 

LD = 3847-5 ” 

WR2 = 8080000-00 (LB-FEET 






in moment of inertia of rotating parts will be benificial. 

3«5 Case 6 Pipe Line with Resistance Dominating: 

The system studied has suction pipe 3847*5 ft. 
long and discharge pipe 23085 ft. long with a head loss 
due to resistance under steady-state of 23.4 ft. and a 
gravity lift of 177,6 ft. As discharge line length is very 
large frictional loss is considerable, and naturally this 
case may bo called as resistance dominating. The equations 
involved and method of solution is exactly same as in case 5a 
or case 5b. 

3.5.1 CasG 6a Small Msment of Inertia of Rotating Parts: 

The graphs plotted for this case are, (h vs t), 

(v vs t), and (a vs t) as shown in Fig. 3. 9, which represent 
transient conditions as a result of power failure. It can 
be observed that the maximum transient head equals the steady- 
state head, 

3.5.2 Case 6b Large Moment of Inertia of Rotating Parts: 

In this case all the data taken are same as for 
case 6a. except the moment of inertia which has been doubled 
to 8080000pound-ft, squared. The method of solution is same 
as for case 6a and equations involved are also tlie same. 

The graphs plotted for this case are, (h vs t), 

(v vs t), and (a vs t), which represent transient conditions 
as result of power failure as shown in Fig. 3,10. It can be 
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CHAPTER 4 


ANALYSIS OF THE TRANSIENTS OF THE TRACY PUJIIPING PLANT WITH 

AIR CHAMBER 


4,1 Introduction: 

An air chamber is a vessel having compressed 
air at its top and connected to the discharge pipe.Fig.4.1. An 
orifice is provided to restrict the inflow into or outflow 
from the chamber. Here a differential orifice has been used 
such that it produces more head loss for inflow into the 
chamber than for a corresponding outflow from the chamber. 

The outflow from the chamber has been made free so that 
column separation in the pipeline is avoided and inflow has 
been restricted to reduce the size of the chamber. 

Analysis of the system with air chamber has been 
carried out because air chamber has its own advantages over 
other auxiliary equipments such as surge tanks. They include 
the following: 

1) The volume of an air chamber required to keep the 

maximum and minimum pr'-ssures within the prescribed 
limits is smaller than that of an equivalent surge 
tank. 

An air chamber can be installed with its axis parallel 
to the ground slope. This reduces the foundation costs 


2) 
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and provides better resistance to both wind and 
earthquake loads. 

3) An air chamber can be provided near the pump, which 
may not be practical in the case of surge tank 
because of excessive height. This reduces the pressure 
rise and the pressure drop in the pipe line. 

4) To prevent freezing in cold climates, it is cheaper to 
heat the liquid in an air chamber than in a surge tank 
because of smaller size and because of proximity of the 
pumphouse. 

The main disadvantage of air chamber is the nocessiv 
to provide air compressors and auxiliary equipment, 
which require constant maintenance. 

A ratio of 2.5il between the orifice head losses f 
the same inflow and outflow has been used. Six cases have 
been studied varying different parameters, which 

CO is initial volume of air inside the chamber, diameter of 
the orifice, water level in the chamber- In the analysis 
of these cases with different parameters, location of air 
chamber has been taken as downstream of the pump, after a 
check valve. 

4.2 Boundary Conditions for Air Chamber; 

Referring to Fig. 4. 2 , at the junction of the 
chamber with the pipeline following equations are availabe: 
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Positive characteristic equation for section( i,n-fl): 


= Cp — Ca^.Hp. 


(Eq.4a) 


i»n+l i,n+l ^ 

Negative characteristic equation for section (i-f-l,l)j 


— Cn -f- • Np 


i+ia i+1,1 

If the loses at the junction are neglected, then 


(Eq.4.2) 


Hp = Hp 

i > n+1 i+1 , 1 


(Eq.4.3) 


Continuity equation: 


i , n+1 


■i+1,1 


■‘■^Porf 


(Eq.4.4i 


in which = flow through the orifice and is 

considered to be positive when flov; is into the chamber . 

Assuming that the air enclosed at the top of 
air chamber follows polystropic relation for a perfect gas, i.e. 


“Pair * ^Pair 


(Eq.4.5) 


in which and v 


T. • and Vr, • are the absolute pressure 
i air Pair 


head and the volume of the enclosed air at the end of time step. 


m is the exponent in the polystropic gas equation. The value 

of constant C is determined from initial stready-state condi- 

m * 

tions, i.e., C= CO . . The values of m are equal to 

1.0 and 1.4 respectively for an isothermal and for an adiabatic; 
expansion or contraction of air. For design calculations, 
an average value of m = 1,2 may be used because transie^it are 
usually rapid at the beginning and they are slow' at the end. 
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The orifice losses may be expressed as 

^^Porf ^orf ^Porf I^Porfl (Eq.4.6) 

in which ^orf = coefficient of orifice loss 

For the enclosed air volume following equations may 
be written: 

“Lir = “p. ^ + W’Porf (Eq.4.7) 

1 p n+1 

^Pair = ^ AC(Zp-Z) (Eqv4.8) 

Zp = Z+O-StQorf+W) f (Eq.4.9) 

in v/hich = barometric pressure head, 

AC = horizontal cross-sectional area of the 
chamber p 

Z and Zp= heights of the liquid surface in the chamber 
above the datum at the beginning and at 
the end of the time stepp 
and CO = volume of air at beginning of time step. 

From Eqs, 4. 1,4. 2 and 4,4 

°Porf = (Cp-Cn)-(Ca.+Ca.^p) Hp (Eq.4.10) 

i , n+1 

and from Eqs. 4, 5,4. 6, 4.7 and 4.8. 

(Hp. ^Hb-Zp-^orfOporf.lOporfl^ [ .C0-AC(Zp-2)]'” =C 

i,n+l 

(Eq.4.11) 
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Eqs. 4,9, 4,10 and 4,11 have three unknowns, 

namely QpQj.fr Hp and Zp. Substituting two of these 

i , n"4“l 

equations one equation in one unknov/n can be derived. For 
example the equation written in terms of, and other 

known variables is as follows. 


Cp-Cn-Qp^^^ ^ 


O.SxD T /p p N 
AC '‘^^orf'^^Porf^ 


~ *^orf ^ *^Porf ^ ^ '^Porf I ^ ^ ^ 

[ CO-0.5 X (Qorf-^Qporf^ ^ - C = 0 (Eq,4.12) 


The corresponding derivative of the function, FU!\[ 
with respect to QpQj,p is given below; 


DFUN 


*i'^'-'*i+l 


^ 9.y Ai^.Q Xj. 2x C X i 0 I 1 X 
^ AC ^ ^orf^i ^Porf • ^ 


[ C0-0,5xDT (Qo^f+Qpo:cf^'" 


Cp-Cn-Qp^^^ 

^i+™^®i+l ‘ 


O.SxmxDT + 1-1^+^ Q.^S1 (Q^^^p+Qp^^p) 


^orf^ ^Porf ^ !^Porfl ^ ^ 


[ C0-0.5xDT (Qorf+Qporf^^ 


( cq ,4,13) 


Then Newton-^Raphson method can bo used to evaluate, 
the value of the unknown Qp^j^f The flow chart for this 

procedure has been shown in Fig, 4,5, 
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4,3 Parametric Analysis of Air Chamber: 

Transient problem studies under case 7. has an 
air chamber just downstream of pump after check valve. The 
length of discharge line, moment of Inertia of rotating 
parts and other datas are same as in case 1. 

It has been assumed that check valve closes 

simultaneously with power failure, so pump characteristics 

vanish from calculations. The boundary conditions, at the 

pump with air chamber has been developed in section 4.2. 

The equations are simplified to have one unknown 

as defined by Eq. 4,12 and derivative of this equation 

has been defined by the Eq, 4,13. Newton-Raphson method 

has been used to find the root of the equation By 

back substituting in Eq, 4.11 we get the value of 

Hp . Zp can be obtained by using evaluated value of 

i , n+1 

Q_, o in Eo, 4.10, For the calculation at next time step 
is assigned the value of Qppppt 4 is assigned the 
value of Zp, and CO is assigned the value of 'Vp^^pp. This 
procedure is repeated upto the required time. This case 
has been studied under six subdivisions by changing the 
values of different parameters involved in the computation. 
The different parameters involved are, the constant 
in which CO is initial volume of air in the chamber, diaineter 
of the orifice and initial water level inside the chamber. 

The input data for different cases are shown in Table 4,1, 
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4.4 Discussion of Results; 

Graphs have been plotted between water surface 
elevation inside air chamber and time as shown in Fig. 4. 2 
and Fig. 4.4, It can be observed that in case 7c and 
7e the fluctuations in water surface elevation are sharp and 
quick. The fluctuations about mean in case If are minimum 
and also the fluctuations are not intense and sharp. 

It is observed from Table 4.1 that the effect of 
2 CO a 

decreasing £ — is that it increased upsurge and 

dov/nsurge. The change in orifice diameter from 7.5 ft. to 
5.0 ft. did not result in any change in the upsurge and the 
downsurge. The effect of decreasing the initial water surface 
elevation is that, maximum water surface elevation decreases 
and minimum water surface elevation also decreases, The 
minimum water surface elevations obtained in case 7c and 
case 7e are negative so chamber designed using parameters 
of these cases are unacceptable. In case 7b minimum cover 
of water column is less than half the diameter of discharge 
line so parameters used in this case are not acceptable. The 
total height of chamber in case 7f. comes to be lov;est 83.4 
ft, but inside air pressure comes to be higher than that 
in case 7d. in which case height of chamber comes to be 

98.4 ft. The total height of the chamber, maximum upsurge 
and minimum downsurge obtained has been given in Table 4.1. 



Case 7a 
















TABLE A*lf Basic Data And Results for different cases with Air Chamber 
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CHAPTER 5 


SUr^RY , CONCLUSIONS AND SUGGESTIONS FOR FUTURE STUDY 
5*1 Summary; 

The Tracy Pumping Plant consists of six large 
units with each pair of pumps connected through a wye branch 

to a 1 5 ft. diameter, reinforced concrete pipe which extends 
approximately 1 mils to the upper canal. Each pump is driven 
by 22,500-hp motor and delivers water at the rate of 767 cu-ft. 
per sec. to the canal under a rated head of 197 ft. On the 
discharge side of each pump there is a butterfly valve which 
closes at a variable rate under the action of servomotor when 
power tripout occurs at the pump. 

A comprehensive study of the hydraulic transients 
for parallel pump system for the Tracy Pumping Plant and for 
different modifications of the original plant involving a total 
of 14 different cases has been carried out. For this a computer 
program has been coded in FORTRAN- IV and run on Dec 1 090 
system at I.I.T, Kanpur. The system has been analysed without 
a valve, with a butterfly valve at the pump, a non-x-eturn valve 
at the pump, a combination of a non-return valve at the pump 
and an intermediate non-return valve at mid section of the 
discharge line. It has also been analysed with some suction 
length, increasing the length of discharge line, and doubling 
the moment of inertia. The system has also been analysed with 
illn airchamber with a check valve upstream of the chamber, by 
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Changing the different parameters involved in the computa- 
tions* All the cases, being analysed under this thesis, for 
the transient problems, are due to sudden power tripout. 

It has been found that maximum head is 1*4 times 
the rated head \'\rfien, system has been analysed without valve. 

The system, analysed vdth butterfly valve has shown that 9' 
percent reducation in maximum head is achieved. The transients 
for non-return valve at the pump are so severe and quick that 
special precautions will have to be taken. In the case of 
resistance dominated system, ilie maximum head rise is equal to 
steady-state rated head, and column separation may occur with 
subsequent rejoining and development of high heads. In case 
of gravity loading, doubling of moment of inertia has decreased 
the maximum head and increased the minimum head, and so increase i 
in moment of inertia will be benificial. The system studied 
with air chamber has indicated that important design parameter 
for the design of airchamber is , the values of upsurge 

and downsurge decrease with increase in this constant. The i 

0 rific o diameter does not seem to have appreciable ef fect on 
upsurge and downsurge. 

The findings of the present study provide valuable 
suggestions for effective and rational design of surge controll-| 
ing equipments. ■ I 
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5.2 Conclusions: 

The following conclusions emerge from the study: 

e) Pressure surges in pump discharge lines subsequent to 
a power failure at the pump motors can be determined 
with satisfactory accuracy if the complete characteristics 
of the pump and control valves are knovm. Only the pump 
characteristics for the zone of normal operation were 
obtainable from the pump-performance data supplied by the 
manufacturer. . This permitted an accurate determination 
of the water-hammer effects in the discharge line up to 
the point at which the flow reversed through the pump. 

However, these data are sufficient for determining the 
possibility of water-column separation in the discharge 
line. Pump characteristics for abnormal operation and 
the characteristics of the butterfly valve are estimated 
for the remainder of the water-hammer solution. 

b) Computation based on the method of characteristics are 
made to study the effect of intermediate non-return valves 
on surge pressures. Effect of valve location, and multiple 
valves are studied. One conclusion of practical interest 
is that there is no strong case for providing intermediate 
non-return valve. 

c) Since the transients are so severe for complete failure 
with check valves for gravity loading, special precautions 


should be taken 



d) In case of gravity loading , doubling of moment of inertia 

decreases the maximum head and increases minimum head ^ 

that 

so it can be concluded ^increase in moment of inertia for 
gravity loading is benificial. 

e) In case of resistance dominated loading maximum head deve- 
loped is equal to steady-state head but the minimum press- 
ure heads show that column-separation may occur with 
subsequent rejoining and development of high heads. 


f) The important design parameter for the air vessel is 
2C0 

, where CO is initial volume of air, if this parameter : 
is increased, the values of upsurge and downsurge decrease, 
also the maximum air pressure decreases and the minimum 
air pressure increases. 


g) The reduction in diameter of the air vessel causes' incre- 
ase in the upsurge and dovynsurge. 

h) The change in the orifice diameter does not seam to have 
any effect on the surges for the small change that has 
been considered in the present study. 

5.3 Suggestions for Future Study: 

The following suggestions are made for further 


study on the basis of the present study: 

a) Effect of a butterfly valve upstream from the air vessel 

with certain specified closure characteristics may be 
analysed. 



Non-return valve v/ith provision of bypass and its size 
effect may be analysed. 

Non-return valve effects with delay in valve closure and 
random delays in the closure of valves may be studied. 

The effect of pressure-regulating valvesj, and pressure- 
relief valves may be studied. 

A comparative study of anairchamber and a surge tank in 
a pumping system may be carried out. 

A program based on the same general concept as for 
parallel pump system may be develop ed for series pump 
system and can be tested for more reliable and correct 


data 
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«a,tv,dxt,tauo,tauf, (¥ ci) 

r 'i'iUMnFR OF POINTS ON TAU VS TIME CURVE =',I2/8X, 

OPLRATxnN ITME = , E8 , 2 , 'SEC ' /8X , 

JN’^’ERVAl for St-ORING TRU curve = ',F6,3 'SEC'/8X, 
TaH VAEUE5'/RX,7F8.3/8X,QF8.3/) 


hEAHlHi, A'>D W«I^.ING INPUT DATA 

IF’’7cASErKurLy**Gn'’To"rO0o"’*"’ 

HEAP C 20 , AHS,ns,bS 


DATA FOR aTRCHAMBER CALCULATIONS 


read 1 . 20 ,*) HSTO, Cn . HB , AM , AC, ZL ,DAC, CORF, QRDI A , CONST .ARORF 
WRi'i’EC3.l, lOliJ USt 5 , CO, HS, AM, AC, ZL,UAC, CORF, ORDIA, const 
r'nhMA'r(8i, 'STEADY STATfc ATRCHAM8ER ABSOLUTE PRESSURE HEAD =',F7,2 
1/8X, 'STEADY-STATE AIR VOLUME ENCLOSED IN AIRCHAMBER a',F12.2 


?/hX, 'BAR nMETRIC PRESSURE HEAD 


,rF5,2_ 


w if ri .1, t ij* • r 

3/BX 'exponent in THE POLYSTROPIC GAS EQUATION s',F5.1 

4/8X,'horizontat, cross-secion al area of the AIRCHAMBER s', F lo. 2 

5/8X 'HEIGHT Op THE LIQTD SURFACE IN THE AIRCHAMBER ABOVE THE 
P datum =' F 7.2 

7/8X, 'DIAMETER OF THE AIRCHAMBER =:',F8,2 

8/8X, 'COEFFICIENT OF ORRIFICE CORF =',E14,2 

9/8X 'DIAMETER OF THE ORRiFiCE =',F15,l 

10/RX,'ThR VALUF of CONSTANT ( 2 . MCOMa/ CQ0*L ) ) =',F5,1/) 


READt20f T) NP,NRLP,IPRINT,NPP,Q0,N0,TLAST,HSUC,TSTRT 
NRlTFC2l,t9) Np NRLP,Qn,NO,TLASf,NPP.HSUC 

FOhMATCRL 'NUMBER OF PIPES =^, I 3/RX , ^NUMBER OF REACHES ON LAST 
lPlPEs>lV8X, 'STEADY STATE DISCH, s' F8.3,'FEET3/S'2X,'STEADY 
2STATE PUMP SPEED s' F6.i,'RPM'/8X,'TIME FOR WHICH TRANS, 

3STATE COND.ARE TO BE COMPUTED s' ,C5.1 , 'SEC'/8X, *NUM8ER,0F 
4PARALLEL PUMPS s' , I3/8X, 'SUCTION SIDE RESERVOIR LEVELs' , F3, i /) 


reading and writing of PUMP data 


HEADC20,)f) NPC,DTH,Qr,HR,NR,ER,WR 2, (FHCX3 »I»lfNPC) 

KEAnC2Q{>J;) (FBCD.isi NPC) 

WRITE(21,23) NpC,DfH,6RrHR,NR,ER,WR2, CFHCI3 ,l*l,NPC) , ^ ^ 
FORMATCSL 'MUMBeA of points on CHARACTERISTIC CURVE s' 14/ 

IBX, 'THETA TNTRVAL FQR STORING CHARACTERISTIC CURVE s', 
2F4.0/8X, 'RATED OISCH, s' F7,2,'FBET3/S'/aX,'RAT£D HEAD s', 
3Fe,l,'FCEI'/8X, 'RATES PUMP SPEED s ' , F6 , 1 , 'RPM'/8X, 'PUMP 
4EFFiSlENCY s' 6j3/8X, 'WR2 s' ,FiO,2, 'LS-^EEr2'//8X, 'POINTS ON 
5UEAD CHARAC'/CX,!oF7,3)) 

W«XTE(21,25} CFB(I)il = li.NPC3 

FORMATC/Sx, 'POINTS ON TORQUE CHARACTERISTIC'/ CX, 10F7,3)) 




data 


FOR PIPES variables 




KEAD{2U,») CL(1),DCI) ,A(I),F(I),IsI,NP) 

21' 4tl5 ' ' ' 

F0RUATC/§X,'PIPE NO' SX.'LENGTH'tBX.'DIA'.BX.'WAVE VEL,' 
J .yX'FRIC FACTOR '/2ix,^CFEEf )'6X, 'FEET', 7X/CFEET/S) ' 5 
wf?iT§C2t,50) ClA^^a),Dh),Aq5,FC|),Isl,§P) 
FOHHATCIoX,I 3,6X,F7,!,3X,F5,2,5X,F7,1,11X,F5,33 
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nT = T,(f'p)/(HR!,»«cf\(f!P) ) 

'•wrTKf 21 

Knj^M,Vrc/§X, 'Plt>F NO', 'iX, 'ADJUSTED WAVE VEL' , 26X , ' CFEET/S) ' ) 
CALCUl.ATinj'oF PIPE CONSTANTS 

pn“fu, 1 = 1, NP ” 

APCT) s:0,78S4*=DCI)’t'*2 

An!vAL\1 = A(l 3 

AW=!l.(T)/(nT*A(I)3 

nCDsAN 

fJUsMfl) 

XF( (A^'-Afjl ) .GF.O.S) NfI)=N(I) + l 
/a-0=li( r)/CDT’»^NCi3) 

WPlFCfkl ,SN) i’,ACT3 
KnHMAT(loX,I3,l2X,F7,U 
CACI)5:G=«‘AR(TJ/Af I) 

C:A’P=CA(I) 

CECT)=:FCl)=^DT/t2,*DCI)*ARf in 
rCll=FC3:3»t,(I)/f2.*G*DCl)*NCI)»ARCI)>K*2) 

C^NTI^FUE 

computation of^cdhstants for pumps 

’THE"FoEr,OwiNG””cONSTANTs”ARE"FOR*Sl"'ufiiTSjFaR"EWGLISH UNITS, 
RFFLAFE 91&04.9P RT 595,875 AND 4,775 BY 153,744_^ 

TR=C595.875+HH»OR)/(nP*ER) 

CS=oaci)»HH 

C6=-a§3.744*TR*DT)/CMR*WR2) 

ALPHAsNO/NR 

V=Qn/(NPP*QR5 

DVsO.O 

OALPHAsO.O ^ ^ 

IF CCASE.DT.5,) GO TO 1006 
C7=NP»CCKci5fCACl) 3*QR/CCA(l)tCACl)*HR) 
Cl5=0,yl3=*'L.S/CARS»ARS*DS>l«8.fG3 

CAEcuLATinN"Dr"sTiADFsTATi'”’cOMDiTlQNS™^^^^^^][[I*.,., 

nORlr-QO 

CinR = QORl + 36,5A 

IF (CASE, EQ, 8, 3 aORl=0, 

IF CCASE,EQ,8.3 QnR=0,5 

C=HST0:<<C0*4>AM 

IF fV,EQ,0.0) GO TO 65 

THs:ATAN2CAUPHA,VJ 

TH = 57,296=1'TH 

GO TO 68 

rH=o,o 

call PARAR(TH,1,Z3 
Hns!Z*HR»(ALPHA»?2 + V**2) 

H(l,i)=HO+HSUC 
HV = H(1,1)-DELH 
call PARABCTH 2,Z3 
BETA=Z«CALPHA»»2+V*=t:2) 

1)0 RO 1 = 1 ,MP 
NOsNCD + l 
DO 70 

H(2,i3«RCl^ft3 


QCl,j3sOO 

QV:36p/2. 


HCI,J5 
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fj'i, 

0 3 

016 

018 

0 


Oil 

015 

Ou 

02 


6 

4 

017 

a 

50 


AOp^an/2, 

40;-, isOGOC 

Af',s7sAHSi-C1 

A!Us3 = 4riS2-Cl 5»Qn=(tOa 

AHi;'i = !\rif>3-c:i.5’F;in^oo 

,V'.S5=4HK4"Ct 5’'Qnf:rja 
COf^TlOijK 

, t, j+nRbH 

Ai*>=Ch(l , D- 1. )/0K 

cn^TtMUE 
aNsO(!''p) + t 
nRCr>sH(Mp,Nf1) 

T-0.0 

WRlTEf^l^HSl 

FnpPA"’c/3X, 'TCMF',2X'ALPHA',4X,'V',4X,'PTPE%17X- 
l'HrAl)(FKSn % /24x/?i!0%5X,'Cl3',6X,'C2)',4X,'(35',5X, 

wprT’p:(24,4oi } 

FOR M at 1 5 X , ' T X Me ' 2 5X , ' DI SCH . FEETi /SEC ) ' , / 1 7 X , ' C I) ' , ox , ' OV % 7X , 

l'C3l ',8X,'(4)',i0X, *C5)'f8Xf 'C63'#8X,'CN + Ur) 

wniTE(22.403) 

FORmaTCSK, 'TIME'. 25X, 'SUCTION HEAD { FEET ) ' , / 1 7X , ' ( AHSl ) ' , 
14X,'CAHS2) ',5X, ^(AHS3)',4X,'{AHS4)',6X, 'CAHS5)*) 
WRlTECil,taif») , ^ , 

FnRMATc/9X. 'TTME',5X,'ZL',13X,'CO',10X,'QOR',6X,'HSTO',/3 
WPiTEf 31 ,i6i«) T,7.L,Cn,Q0Rl,HSTn 
FnRMATC/SX,F7.1 ,3X,F7.2,3X,F12,2,3X,2F9,2) 

K=:o 
1 = 1 

NNaNCT)+l 

IF (CASE, ED. 4.) GO TO lOll 

WRITEC21,8&5 T AGPHA,¥,1,HC1,1),HV,HCI,2),HC1,33,HC1,4J, 

IHCI ,5) ,H(1 ,Nfl) ,AD 
IF CCASE.NE.4,) go to lOlS 

WRITE (2 1.843 T, alpha, V, I, HC 1,13 ,H ( 1 , 2) , H C 1 , 4) , H C 1 . 6 3 , H ( 2 , 1 ) , 

. AH , . , 

WRITE(22,4023 T . AHSl , AHS2 , AHS3 , AHS4 , AHS5 
FnRMATC2X|F7.1 ,2x,5Fil .23 
WRlTEOl.iOl?) T,ZL,CO,OOR,HSfO 
00 B9 1=2. HP 
NM=N(T)+1 

F0HMAT13F7,2.2X,I2,2X,7F8,1,2X,F8.23 

FnRMATtr7,i,2F7.2.2X,t2,2X,7F8U,2X.F8,23 

FnKWAT(5X,F7,l,3X.F7,2.3X,Fl2,2.3X,2F9,2> 

continue 

T»X+DT 
K=K + i 

IF CT,GT,TLAST3 go TO 240 



IF CCASF.EO,?.) ALPHA* 0,0 
IF (CASE, EQ, 7. 3 V= 0-0 

CN = 0 ( 1 . 23 -Ht 1 , 2 )*CAC 1 3 -CFCl 3 *Q(i. 2 )*ABSCQCl , 23 ) 
IF (CASE. LT, 7, 5 GO TO 1 H 2 
CALF GUIDE(ANS) 

IF CANS, LE. 0,03 CORFxCORF /2 ' 

CALF CUIDKCANS) 

IP (CASE. EQ. 8,3 CALL PUMP 
IF CCCASE,eO, 8 ,),AND,CV,LE, 0 |O 3 ) V= 0,0 
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1 !=' ),ANn, fv.LR.o.O)) alphasO.o 

Of’C 1 , 1 )=«A.'Jf. + MpPi»^VtOR 

7,P = 7,LH 0, C OOH+ Af4f> ) *0'^/AC 
VP=rn-Ar*czP"?.L) 

ZL=:7F 
CP = ^'F 

IF fC'VSF.FQ.e. ) GO TO J1J3 
IF fCAse.EO,?.) GH TO 1113 
CAIJ. PiJ.’^P 

IF CCFASK.tO, J.I.AND.CV.LF.O.O)) V=0,0 
IF f CCAffO.feiP.i. ■), AND. CV.EQ. 0,0)1 ALPHAsO.O 
IF f fF.vSE.E0.A.).ANO,(V.LE,0.0)) V=0.0 

IF rcrA.Sfc:.c:n,4,).AND,(v.EO,o.O)) alpha=o.o 

>.«|PC S 1 )"i<PP^VFQR 
uvznpf 1,1)/?, 

A n p S!: 0 V 

CR=rw/C5 

HPci )-cop(i,i)-rN)/rAcn 
HFnPsrf)PF*AHr»’^ABS('AMS) 

HPrtWsHPCl ,l)+-HB"ZL-HpnR 
HBTHsHPAr 

IF CT.GE.TSTRT) HPCI ,n = COP(l,l)-CN)/CACl)-CV=^QPCl,l) 

l»ARSfOP( 1 , 1) )-8, 17 

IF fCASE.EQ.l.) CV=0.0 

IF CCAse,EQ.3.) CVsO.O 

IF CCASE,Ea,4.) CVsO.O 

IF (CASE.FO.S.) CVsO.O 

IF CCASR,EQ,6,) CVsO.O 

IF CCASE,FQ.7.) CVsO.O 

IF CCASE.FQ.S,) CVsO.O 

HV=HpCl,i)-CV*QPCl,l)*ABSCQPCt ,1)5 

Ah=fHP(l,l)-4.)/HP 

IF CCASE.lt. 5.) GO TO 1003 

IF CCASE.FQ,?.) go TO 1003 

AHSlsHSOC 

AH52 = AHS1-C15*QP(1,1 )»ABSCQPCl,n ) 



on 170 1=1, NP 

NflsNCl) 

on 160 0=2, HN 

CM=ACI,J + 1 )-CACI5>*'H(I,J + l)-CF(l5*QCI,J+l)*ABS(QCIfa+13 
CPsOCT, J»l)+CACl)*HCl,J-l)-CF(I)’t'QCI,0-l)*ABS(OCI, J-15 
aPf.i,a)=o. 5 »cCp+cN) 

HP(I, J)aCGP«QPtT,a))/CACI)-CV*QPCI,J)*ABSCQP(T,a5 ) 
IF^(|^|E,TSTRT) HPCT,J)=(CP-QP(I,J))/CA(I) 

'ConTlNHE 



ilPlssNP-l 

iFCNP.EO.t) GO TO 178 
on 175 1=1, Mpi 
MlsPCI^ 

Nli*N(I)4-l 


n t w If i A j V 1 

CNaOaTl, 2 )-CA(I + l)»HCI+i, 2 )»CFCI + l)*aCI + l, 2 )='‘A 8 SCQ(lTi, 23 ) 


Cf=0CT,^'l )+^ACl)=«'HCI,NU-CFCI)*QCT,NntABSCOCI,Nn) 

!i»(I ,M».;)sCCP«f'N)/CCA( T)+CACI + U) 
t + U ,WW) 

IF ((CASE.L-'J.^, ■J.Af^O.CQPC r,NN).E;Q, 0 , 0 )) HP( I , NN 5 aCP/CA ( I) 

IF CCCASi:;.f/"j, 4 .->.AND,CQP(T + l,l),E 5 . 5 . 0 )) HP C I + 1 , 1 ) a-CN/CA C 1 + J ) 

caM-uN'is 

REoEftVUTR AT UOWNSTREAH "eMD : 

W'vsNpiP) t.l 
HPtMP,NN)=HRES 

CP= 0 C«P,NN -1 ■j+cACHP)^=H(NP,NN« 1 )-CFCHP) 4 :QCNP,NM- 1 H 
i ABSC 0 (WP,MW- 1 ) ) 

QPC*^ 5 P,HM)=CP'-CaCMP)*HPCNP,NN) : 

CParp/c 5 

vSTQRiNG MAX, AND MiN, PRESSURES AND VARIABLES FOR NEXT TIME STEl 
DO 730 1=1, wp 
NNsNCI)+l 
on 22 n J= 1 ,UN 
aCl,J)=:QpCl , J) 

HCI,J)=HP(I , J) 

IF (HCI,! l.GT.HMAXCin 
IF (HCl / n .t.T.HMlNCI) ) 

CONTINUE 

CONTINUE 

IF CR.EQ.TPRlNf) CD TO PO 
on TO 150 
WRITE(21,250) 

' .PRESS. ', 5 X, 'MIN. PRESS.'/ 


HMAXCI>=H(ia) 

HMIN(l|=H(I,l) 


FOrMATCZ/IOX, 'PIPF N 0 ,', 5 X,'MAX 
12 SX '(FEET) % lOX, '(FEET) '/) 

WRITE ( 21 , 260 ) (f,HMAXCl),HMIN(I),T=l,Np) 
FORMAT( 12 X,I 3 , 7 X,Eni, 9 X,F 7 ,l) 


STOP 
END 

subroutine 


PUMP 


dimension FH(80) ,FB(80),y(S0) 

common /CP/ alpha, QR,V,CN,DALPHA,DV, BETA, C5,C6,NPP,T, CASE 
COMMON /PAR/ fh,fb,oth,tauo,tstrt,hr,qo,ho,cv,hsuc 
COMMON /RAH/DXT,y ,C7,C8 
KKrO 
JJ = 0 

’"’"'computation"”"*!? pump " discharge 

VE=V+DV 

ALpHAEaALPHA+DALPHA 
J0=JJ+1 

IF CVE.EQ.O.O.AMD.ALPHAE.EQ.O.O) go to 29 
TH=ATAN2(ALPHAE, VE) 

THi=TH 

TH s TH»57.296 
IF CTH.LT.O.O) THaTH+360, 

IF CTHl.Lt.O.O) TH1STH1+6.2831S 
GO TO 30 
TH«D,0 
TH1=0,0„ ^ 

MaTH/OTH+1 

Al=FH(H)#M-FiI(M-H)’»‘(M-l) 



86 


A?=:CFT»C»'' + 1 )/(DTH*0,0174533 

A3sPrtfM)?'’-.FeCM + 13 ♦CM-l ) 

A4sfFP(M+1 )-FBCM) )/cdtH’<' 0. 0174535 
Al.'PSQsALpHAFf alphas 

i/FyOsVe^'VF 

aLpv=alpso+vesu 

call PARftBl fT,TAiJ) 

IF CCASF.Ey.^.) Sn TO 1004 

IF fCASF,R.3.l.) Sn TO 1002 

IF (CASE, SO, 3.) GO TO 1002 

ir CCASF.Ka,4.) GO TO 1002 

IF (CASS, Fa. 5. 3 on to 1002 

IF ccASF.ca.e, ) go to 1002 

IF (CASE, Fa. 7.) GO TO 1008 

CV = C!.0 

Cn = (CAl*C54CP’*'CS5/(C5*C53 

CVsB.l 7*2.»VE't'VF*WPP»fJpP*HR/CTAO>»=TAU^OO»QO=»=C5) 
Fl.=r5’»:A14ALPV + C5=l'A2*ALpV*THl-QR«yE*NPP+CN4-CS»HStJC/HR 
IF (CASF.FQ,7.) GO TO iOtJP 
IF (CASS, FQ. 1,5 F1=F1-C5»HS0C/HR 
IF CT.GE.TStRt) 

1 Fl=Fl“C"i*Ctf*QR*OR*VF*ABSCVE)*NPPI'MpP/HR 
F7=ALPHAE-C6tA3^ALpV-r6*A4=t^ALPV*THl-ALPHA-C6=»'BETA 
Ft AL=C5*C2.»Al»ALPriAE+A2»VE+2,»A2*ALPHAE>t=THl ) 

FI V=:C5#(2.*A1»VE-A2»ALPHaE+2.»A2»VE*TH13«QR*NPP 
IF (CASE.EQ,?.) go to 1010 
IF CT.GE.TSTRT) 

1 FlV«FlV>-2»C5»C’V=('QR*QR=«'NPPfNPP’t‘ARSCVE)/HR 

IF (CASE.EQ, 5. 5 Fl =Ai ♦ALPV+A2*ALPV’t:THl-C7»VE+C8 

IF (CASE.EQ. 6,) Wl^M*Ahpy+k2*hhPV*THi^C.l*VE+CB 

IF (CASF.EQ.5.) FI 7=2 , ♦ A1 tVE+a , ♦A2*VE*TH1-A2* ALPHAE-C7 

IF (CASE.EQ. 6.) FI V = 2,=»'Al+VE+2.>t'A2*VE*THl-A2*ALPHAE-C7 

IF (CASE.EQ, 5.) FlALs2.»Al=<‘ALPHAE+2,*A2*ALPHAE*THl+VE*A2 

IF (CASE.EQ. 6, 5 FI AL=2 , ^Al » ALPHAE42 , * A2*ALPHAE»TH1 ♦VE't'Aa 

F2aL=!1.-C6*C2,*A3*ALPHAE+A4*VE+2.*aUaLPHAE»IH15 

F2V=Cfi*C-2,«A3*VE + A4»ALPHAE-2.#A4*VE4<THl) 

DENnM=FlAL»F2V-FiVi»!F2AL 

DALPHA=CF2*F1 V-F1*F2V)/DEN0M 

DV = (F1»F2AL-F2»F1AL)/T>EN0M 

AT,PHAE=ALPHAE+DALPHA 

VE=VE4DV 

IF (Aas(DV) ,bE. 0.001 .OR. ABSCDAIiPHA3.LE. 0,001) GO TO SO, 

IF (JO.GT.iO) GO TO 70 
Go TO 8 

TH=ATAM(ALPHAE/VE) 

TH=ATAN2(ALPriAE,7E3 

TH=57,296#rH 

.(F (TH.LT.U.O) TH=TH + 360, 

CALL PARAR(TH.2f BETA) 

IF (MB.EQ.M) GO TO 60 
HBsTH/1)TH + 1 
(F (.HR.EQ.M) Go to 60 
GO TO 8 

QALPHAsALPHAE-ALPHA 

i)V=VE-V 

nLpHA=ALpHAE 

7®VE 

beta =BETA*CaLPHA»ALPHA+V’*'V) 
return 

Fof<MAT(li!'*»J'fTERATfois IN PUMP SUBROUTINE FAILED' /8X, 
l'T=%F8.?/8X,'ALPHAE =',Ft§,4/8X, 'VE a',F6,33 
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P^kABOLA THRRHuGH three points on curve 
(I) IS the T.VTEOER part OE the ERAC^ION 
(R) IS THE heal part OE THE FRACTION 

GFNFKALf'Y VAf-uES OF fY) ARE KNOWN FOR EOAL INCREMENT IN (X) 
if' VALUES UE (Y> ARE DESIRED FOR (X) WITH ORIGIN (0,0) 

Oi'E FIRST FINDS A PARABOLA THROUGH THE THREE POINTS WITH 
A'-'J-S.. vertical by TRANSFERING THE ORIGIN TO (X',Y'3 
THERFORE Y' s: A*X'’’ 4==<'2 + B*X'’ 

LET XfN + n-X(M):;X(N)-X(N-l)sH 
THEN ry(Nti 3 -Y(N)rA»H»=)' 2 +B>(‘H 3 
hND rj'lN-’J.)-YCN)sA*H*f 2 -B»H 3 
A=[Y(.Nt-]i 3 + Y(u-l)-2.TYf N)3 /C2.*H»*?3 
Hslycn + I )-y(N-l)3 /C2,*H’I'!^2) 

therefore, tY=tYCN) + 6,5*THETA*=»;2CCYCN + 13+Y(N-l)-2,»YCN3) 

+ TI!E'}'A*C (YCN + 1 )-YCN-l))»0,51 3 
IN WHICH THETa^H HAS BEEN SUBSTITUTED FOR CX') 


COMMON /PAR/ FH,rR,DX,TAUO,TSTRT,HR,Qn,HO,CV,HSUC 
DTHEfJS| 0 « FHCBO) ,FBC 86 ) 

IsX/DX 

R 3 (X-T*DX)/nx 
IF CI.EO.O)R!sP» 1 . 

1 = 1 + 1 

IF CI,LT, 2 ) 1=2 
(50 TO ( 10 , 2 R 9 ),J 

0 ZsFH(T) + 0,5*R*CFH(I + l)-rHCl-l)4-R*CFH(T + l)+FHCl-l)-2,*FHCl)) 3 

RETURN 

9R Z=FBCI) + .5«iR»CFPCT + l )-FBCI-t)+R»(FBCI + l>+FBCI-l)-2,»FB(I)33 

RETURN 
F.ND : 

SUBROUtTnE PARABI (r,TAU3 

COMMON /RAM/DXT. Y,C7,C8 
dimension YC603 


I=T/DXT 

rsCT-T^-'DxTJ/DXT 
IF Ci.EO.O.n) H=R-1, 

IF'^CI IjT I?) 

tansy? i);b.f^PtCY(I + U-YCl'-U+R*Cy CI + l)+YCI-13-2.*y(I))3 

r^ETURN 


the FOLLDWING program calls THE newton-raphson subroutinep 




SUbPOUTIMF GL!Tl)E(ANS) 

COMMON /RAMA/ HSTD . CO , HR , AM , C .CORE, AC, CAT , DT , ZL, QOR ^ 
COMMON /CP/ alpha, OR, V,CN,0ALPHA,DV, BETA, C5,C6,!IPP,T, CASE 

external EUN.DFUN 
XO 1=0 0 

cal?, root C fun, DEUN.XOI, 0, 001, 50,10, AMS3 

RETURN 

END 
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5!!hP(.!!Ti?jr-: foh Fj,NDTNf: thf root of m equation ri 

«F^rrir.i-RA,PHSD^' ME'f’HOU. FUN IS THE FUNCTION NAME, 

XT„Tt, IS the XMJTIAL trial ROOT. 

D>"u'' JS •’’fiF DrPJVATIVh OF THE FUNCTON FUN, 

K^S ,IS Tut: FRFsrKTBED ERROR LiHMIT OF THE ROOT, 
rfAXTV IS THF HAXiMUMOM NUMBER OF ITERATIONS ALLOWED, 

INDTC: IS Au IMuTCAIOR SF,T to 1 IF THERE IS NO 
COi'-iVEFGEUCE IN MAXTT ITERATIONS, 

THE ROOT IS RITURNED IN SOLN, 

vSnnRlIilTTNE RUOTfFUNrN,uFUNCN,XlNTL,EPS,MAXIT,INDiC,SOLN} 

cnuMnN /HAMA/ HSTO,CO,HB, AM,C,COKF,AC.CAT,DT,ZL,aOR 
common /C'V alpha, UR, V,C fnDALPHA,DV,BETA^C5^C6 ^NPP,T, case 
iUOlCsO 

DO 100 IsI.MhXIT 
FOsFUNCN(CN,XtNTL,CAT-HR,ZL 
DFo = r)rUNCMCCNrXTNTL,CAT,HB, 

XNKW=XINTL-(F 0 /DF 0 ) 

IF CAHSC (XNE?i-XlNTL)/XNEW) ,LE. EPS) GO TO 120 
XrNTLsXNEW 
U CONTINUE 

i INr;IC=l 

SOLNsXNEw 

return 

0 SOLNsXNEW 

return 
end 



the OFFINITION of the functions ARE given BELOW 
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function FUN(CN,X,CAT.HB,ZL,Q0R,DT.AC,C0RF,C0,A«,C) 

f!!N=CC-CN-X)/CAT+HB-ZL-0.5F(X+Q0R)^DT/AC-C0RF»X*ABSCX)) 

l*(r(l-DT»Q, 5 l:CQnR'lX))**AM-C 

RETURN 

CUD 


FUrJU'CTON DFUN C CN , X , CAT , HB . ZL , QOR , OT , AC.CORF , CO , AM) , ^ 

!)FuM=:-.((-CN-X)/CAT+H8-ZL-0,5*CXTQOR5l'DT/AC-CaRF*X*ABS(X) ) 
l*AMfDT»0.5=t'CCCO-DT=ro,5*COORTX)))#l'CAH-.l) 

2 + C-l/CAT-0,5»DT/AC-"2.0*CnRF*ABS(X))»{CCO-DT*O,5*CQOR+X3)5 


return 

£NU 


THE end 




